Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN MfeK style file vl.4) 



Cluster mass reconstruction from weak gravitational 
lensing 

Gillian Wilson, Shaim Cole and Carlos S. Frenk 

Department of Physics, University of Durham, Science Laboratories, South Rd, Durham DH1 3LE 
. . . 5 February 2008 

\o : 
o\ ■ 

0\ • ABSTRACT 

£h ' Kaiser & Squires have proposed a technique for mapping the dark matter in galaxy 

clusters using the coherent weak distortion of background galaxy images caused by 
gravitational lensing. We investigate the effectiveness of this technique under con- 
CTN ' trolled conditions by creating simulated CCD frames containing galaxies lensed by a 

, model cluster, measuring the resulting galaxy shapes, and comparing the reconstructed 

mass distribution with the original. Typically, the reconstructed surface density is di- 
minished in magnitude when compared to the original. The main cause of this reduced 
signal is the blurring of galaxy images by atmospheric seeing, but the overall factor by 
which the reconstructed surface density is reduced depends also on the signal-to-noise 
ratio in the CCD frame and on both the sizes of galaxy images and the magnitude 
. limit of the sample that is analysed. We propose a method for estimating a multiplica- 

tive compensation factor, /, directly from a CCD frame which can then be used to 
\Q , correct the surface density estimates given by the Kaiser & Squires formalism. We test 

ON ' our technique using a lensing cluster drawn from a cosmological N-body simulation 

with a variety of realistic background galaxy populations and observing conditions. We 
find that typically the compensation factor is appreciable, 1.4 <; / <; 2.2, and varies 
I 1 considerably depending on the observing conditions and sample selection. We demon- 

strate that in all cases our method yields a compensation factor which when used to 
^3 i correct the surface density estimates produces values that are in good agreement with 

those of the original cluster. Thus weak lensing observations when calibrated using this 
method yield not only accurate maps of the cluster morphology but also quantitative 
^ • estimates of the cluster mass distribution. 
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1 INTRODUCTION 

In the standard picture of hierarchical structure formation, 
clusters are the most recently formed bound structures and, 
of all objects in the Universe, they are expected to retain 
most traces of the initial conditions which determined their 
formation. It has become apparent in recent years that the 
dynamical state of clusters can be probed effectively by 
analysing the distortions in the images of background galax- 
ies gravitationally lensed by the cluster potential. When 
suitably analysed, these distortions provide a direct mea- 
sure of the cluster mass as well as a map of the distribution 
of dark matter within the cluster. 

Traditionally, estimates of cluster masses have been 
based either on the virial theorem, or on the properties of 
the hot X-ray emitting intracluster gas or on a combina- 
tion of both (e.g. Hughes 1989). In all cases, a number of 
assumptions are required which introduce unavoidable un- 



certainties. For example, only the radial component of the 
velocity of cluster galaxies is measurable, so assumptions 
need to be made concerning the missing information about 
the galaxies' orbits in three dimensions. In addition, all op- 
tical observations are confused by chance alignments of field 
or group galaxies physically unrelated to the cluster (Frenk 
et al. 1990). X-ray observations are less influenced by pro- 
jection effects, since the bremsstrahlung radiation from the 
intracluster gas is proportional to the square of the gas den- 
sity (e.g. Fabian 1988). However, the inferred mass depends 
on the temperature profile of the gas and this is still poorly 
constrained by existing X-ray data (e.g. Arnaud 1995). Fur- 
thermore, since gas density falls off rapidly with distance 
from the centre, other techniques are required to measure 
mass at large distances from the cluster centre. 

The use of clusters as cosmological tools is not restricted 
to the information provided by their total mass. If small 
systems of galaxies have recently merged to produce a rich 
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cluster, evidence of vestigial substructure should be appar- 
ent. The distribution of mass within clusters is therefore as 
important a diagnostic of the cluster formation process as 
is their total mass. For example, Evrard et al. (1994) and 
Mohr et al. (1995) have used simulations of cluster gas and 
dark matter to suggest that an X-ray morphology-cosmology 
relationship exists. They find that clusters formed in low 
S7 models are more regular and spherically symmetric than 
clusters formed in the = 1 case. This is a reflection of the 
fact that clusters form earlier in low Q universes (Lacey & 
Cole 1993). 

Uniquely amongst all techniques for studying galaxy 
clusters, gravitational lensing is directly sensitive to the dark 
matter within the cluster. Thus, lensing studies bypass the 
uncertain connection between the luminous material in clus- 
ters and the dynamically dominant dark matter component. 
In principle, gravitational lensing provides the most power- 
ful tool available to extract cosmological information from 
clusters. 

The first detections of gravitational lensing by clusters 
were made in the late 1980s. Lynds & Petrosian (1986) 
reported the discovery of giant arcs in the clusters A370, 
A2218 and C12244-02 and, independently, Soucail et al. 
(1987) discovered arcs in A370. Giant arcs are spectacu- 
lar but rare occurrences. Their existence depends upon the 
serendipitous alignment along the line-of-sight of a back- 
ground galaxy with a dense cluster core. Perfect alignment 
behind a spherically symmetric core region would lead to a 
perfectly circular image - a so-called Einstein ring. In prac- 
tice only portions of the ring are produced, causing the im- 
ages to be called arcs. A comprehensive review of giant arcs 
and their properties may be found in Fort & Mellier (1994). 

It is far more common for a galaxy lying behind and to 
the side of a cluster to be stretched or sheared tangentially 
only slightly. Galaxies which have undergone only weak dis- 
tortion are generally referred to as arclets. These galaxies are 
too faint for spectroscopy so individually they are imprac- 
tical as indicators of lensing. However, the cluster mass dis- 
tribution can be recovered statistically by analysing collec- 
tively these weakly, but coherently lensed arclets. Tyson et 
al. (1990) were the first to study weakly lensed images in the 
clusters A1689 and C11409+52. In their pioneering study, 
they observed an excess of tangentially aligned galaxies and 
set constraints on the cluster potential from this data. Subse- 
quently, a number of authors (e.g. Kochanek 1990; Miralda- 
Escude 1991) have attempted to determine cluster param- 
eters, such as velocity dispersion and core radii, from ob- 
servations of weakly lensed galaxies by model fitting. They 
assume a priori some form for the distribution of mass in 
the lens and then determine the most likely values of the 
model parameters. 

Kaiser & Squires (1993), hereafter KS, proposed an ele- 
gant, model-independent mass reconstruction method. This 
technique, described in detail in Section 2.2, produces a 
"map" of the surface density at each point in the cluster. 
Since the initial idea was proposed, progress on the theo- 
retical front has been extremely rapid. C. Seitz and Schnei- 
der (1995a) and Kaiser (1995) have developed extensions 
of the method capable of simultaneously reconstructing the 
cluster mass in the weak and strong lensing regimes. The 
KS technique assumes that lensing information is available 
over an infinite field of view. In practice, the limited size of 



CCD frames introduces spurious boundary effects. Schnei- 
der (1995), Kaiser et al. (1994b) and S. Seitz & Schneider 
(1995b) have all addressed this problem. Some attempts 
have been made to investigate the KS method using clusters 
grown in N-body simulations as lenses. Bartelmann (1995) 
used a sample of 60 clusters from the simulations of Bartel- 
mann & Weiss (1994), synthetically lensed them, and inves- 
tigated a variety of reconstruction algorithms based on the 
KS method. 

It has long been realised that the original KS technique 
does not furnish the absolute value of the mass surface den- 
sity because a uniform screen of mass located between the 
observer and the background galaxies produces no gravita- 
tional distortion. Bartelmann & Narayan (1995) and Broad- 
hurst, Taylor & Peacock (1995) have suggested methods to 
break this degeneracy by utilising the magnification of the 
lensed background field galaxies rather than the distortion 
of their images to try to constrain the absolute value of sur- 
face density. Whether these theoretical ideas can be applied 
effectively in practice still remains to be seen. 

With the widespread availability of large CCDs, obser- 
vational studies of weak lensing have proliferated in the past 
couple of years. Bonnet, Mellier & Fort (1994) detected a 
lensing signal in C10024+1654; Fahlman et al. (1994) and 
Kaiser et al. (1994a,b) in msl224+007, A2218 and A1689; 
Smail et al. (1994) in C11455+22 and C10016+16; and Tyson 
& Fischer (1995), again in A1689. These detections have gen- 
erated a great deal of interest in lensing studies as well as 
some controversy. For example, the dark matter mass in- 
ferred for msl224+007 from lensing data by Fahlman et 
al. (1994) is three times larger than the virial mass in- 
ferred from optical data. (See also Carlberg, Yee & Ellingson 
1994). A2163 is another paradox: it is the most X-ray lumi- 
nous cluster known but no lensing signal has been detected 
(Squires 1994). 

It is clear that weak gravitational lensing is a power- 
ful and useful technique but that there is still much work 
to be done in order to understand the problems and sys- 
tematic effects implicit in its use. The values of the cluster 
mass surface density inferred from distorted images of back- 
ground galaxies are complicated by the influence of inherent 
observational effects such as noise, seeing, crowding and pix- 
elation (i.e. the discrete sampling of the average intensity in 
the detector pixels). The quality of the signal will depend 
additionally on intrinsic properties of the lensed galaxies - 
their magnitudes, sizes, ellipticities and redshifts. The rel- 
ative importance of all these factors will be investigated in 
this paper. 

Our aim is, firstly, to confirm that the KS technique 
does indeed recover accurately a complex lensing mass dis- 
tribution under controlled conditions. We do this by simu- 
lating CCD frames of galaxies which have undergone lens- 
ing and analysing these frames with the same techniques 
of faint galaxy data reduction that are commonly applied 
to real data. We construct a variety of artificial clusters to 
investigate the relative importance of observational effects 
by varying them individually. Since the mass distribution of 
the artificial clusters is, of course, known, the accuracy of 
the mass reconstruction obtained under differing observing 
conditions can be assessed. In general, we find that the re- 
covered mass surface density is less than the true surface 
density by some factor. It is impossible to tabulate a com- 
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pensation factor for all possible combinations of variables 
but by means of an example we illustrate a method for esti- 
mating this factor for any given data set. Several groups have 
applied a correction for the effects of seeing on their data. 
Fahlman et al. (1994) and Tyson & Fischer (1995) did this 
by modelling the properties of background galaxies. More 
recently, Squires et al. (1995) used HST images which they 
artificially sheared and degraded to simulate ground based 
observations. The method described here complements these 
techniques as it does not require any assumptions about the 
background galaxies or additional data. Our intention is not 
to model conditions with any one particular telescope or 
observational set-up in mind, but rather to produce results 
which can be applied generally. 

In Section 2 we summarise the lensing concepts and 
equations which we will require. In Section 3 we describe 
the details of the source galaxies, cluster lens and analy- 
sis software which we use. In Section 4 we generate simple 
spherically symmetric lenses of varying mass in order to il- 
lustrate the effect on the reconstructed surface density of 
atmospheric seeing and of non-linear terms ignored in the 
KS method. We then proceed to show how the reduction in 
the lensing signal caused by atmospheric seeing can be com- 
pensated for by performing a calibration exercise in which 
a known shear is applied to a CCD image. Here we use a 
cluster drawn from a cosmological N-body simulation as a 
lens and perform the complete analysis using realistic dis- 
tributions of ellipticity, size and redshift for the background 
galaxies. In Section 5 we repeat the analysis varying each 
of these distributions, demonstrating the versatility and ac- 
curacy of our calibration technique and the general power 
of the KS reconstruction method. We conclude in Section 6 
with a summary of our main results. 



2 WEAK LENSING 

2.1 Basic Concepts of Lensing 

The basic lensing geometry is shown in Fig. 1. The primed 
letters are points in the source plane. Their unprimed coun- 
terparts are corresponding points in the image plane. Light 
from the source galaxy S' follows the path S'lO to the ob- 
server. The galaxy's apparent position in the source plane is 
I'. As can be seen from the figure, for a circularly symmetric 
potential, galaxies are displaced radially outward from the 
centre of the cluster ie from S' to I'. Implicit in this scheme 
is the assumption that all deflection takes place at one point 
(the lens) in a light ray's journey to the observer. We use 
the subscript s to refer to the true, unlensed position of the 
source and the subscript i to refer to the apparent position of 
the image. The letter D denotes angular diameter distances. 
Its subscripts denote observer, lens and source galaxy, a is 
the angle through which light rays are deflected at the lens 
and (3 is the apparent angle of deflection at the observer's 
position. 

In the small angle approximation these angles are re- 
lated by 



D\ B a = ST = D s [3 



(2.1) 



Hence, the lens equation, relating the true position of a 
source galaxy 9 B to its apparent position 9\ by means of 
the bending angle /?, is given by 
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Figure 1. Schematic lensing diagram 
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The angle, a, through which rays are deflected can be calcu- 
lated from the gradient of the lens' projected gravitational 
potential (see e.g. Schneider et al. 1992, Chapter 5). Thus, 
if we define a dimensionless 2-D potential by 



then 
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I dz 



(2.3) 



P(0i) = V$2D (2.4) 

where V = ( , ) is the gradient operator in angular 
coordinates on the sky. Just as the three-dimensional poten- 
tial is related to the density via Poisson's equation, so the 
two-dimensional potential is related to the projected surface 
density, S, via the two-dimensional Poisson equation: 

V 2 $ 2D = (2.5) 



where S CI it , the critical surface density, is defined as 



$crit — 



4ttG DoiA, 



(2.6) 



This critical surface density is that required to form multiple 
images of a source object. It is also the mean surface density 
within the Einstein ring. 



2.2 The KS Inversion Technique 

A point mass produces a distinctive distortion signal. The 
images of surrounding background galaxies are elongated 
in the direction of tangents to concentric circles centred on 
the point. For a complex lens, if one chooses a point 9 in 
the lens plane, then the correlation of the actual pattern of 
galaxy orientations with the tangential pattern which would 
be produced by a point mass at 9 is a measure of the surface 
density of the lens at that point. Remarkably, as shown by 
Kaiser & Squires (1993), in the weak lensing regime, where 
S <C Scrft j the lens surface density is simply proportional to 
this degree of correlation. Define 
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a(9) = 



s-s 

Sctit 



(2.7) 



where o(9) is the estimated deviation of the lens surface den- 
sity at position 9 from the mean surface density S within the 
area being considered ("over-density" or "under-density" ) , 
measured in units of the critical surface density S^it • Kaiser 
& Squires were able to express the lens surface over-density, 
a, at any position 9 in terms of a direct sum over the galax- 
ies, positioned at 9 g . They showed that 



(2.8) 



The sum in equation (2.8) is performed over all the 
galaxies in the image and n is the mean number density 
of galaxies per unit area. The t hree remaining factors on 
the right hand side of equation (2.8) are the ellipticities of 
the galaxies, d, the kernel, x»; which is the distortion pat- 
tern produced by a point mass, and a weighting function, 
W(9), which produces a smoothed estimate of the lens sur- 
face over-density. Both the ellipticities, d, and kernel, Xi are 
two component quantities and summation over i is implicit. 
We now define these quantities more precisely. 

The components, e iy of the galaxy ellipticity are the 
following combinations of the intensity- weighted second mo- 
ments of the image: 



ei = 



and 



ei = 



2Ix 



/xx + Iyy 

For example, the intensity moment I xy 

_ J F(9)(9 x -9 c x )(9 y -9 c y )d 2 9 
Xy / F{9)d?9 



(2.9) 



(2.10) 



(2.11) 



where F(9) is the intensity at 9 and 9 C is the centroid of the 
galaxy image. For a galaxy whose isophotes are concentric 
aligned ellipses with axial ratio b/a, the size of the elliptic- 
ity, e = yje\ +e 2 2 , is simply e = (1 - (b/a) 2 )/{l + (b/a) 2 ). 
The components ei and ei are given by ei = ecos(2^) and 
62 = esin(20), where <f) is the angle between the x-axis and 
the major axis of the ellipse. Note that for circular galaxies 
(b/a = 1) e = 0, while for highly elliptical galaxies (b/a <C 1) 
e -» 1. 

The components of the kernel x*($g — are given by 
the expressions 



smoothing is assumed then the variance in the KS estimator 
is formally infinite. The smoothing function is required to 
be a low-pass filter but is otherwise arbitrary. The choice 
of window function will, in general, depend on the specific 
property of the lens mass distribution one is interested in. 
Here we wish to make maps of the mass distribution and so 
we have a adopted a simple Gaussian window function with 
transform 



T(k) = exp ■ 



(2.14) 



The smaller the smoothing angle, # sm , the larger is the 
typical error in a due to the intrinsic ellipticities of the 
background galaxies. For this choice of smoothing window, 
Kaiser & Squires (1993) compute the variance in the esti- 
mator, <r, 



((a -(a)) 2 ) 



(2.15) 



where e is the mean value of the intrinsic ellipticities of the 
galaxies used in the reconstruction. The uncertainty in a 
is proportional to the rms galaxy ellipticity and inversely 
proportional to the smoothing angle and the root number 
of galaxies per unit area. In practice, however, uncertainties 
also arise due to errors in measurement, pixelation, noise 
etc. 

We have chosen the smoothing angle 9 sm to be 0.25 ar- 
cminutes as a compromise between producing a high reso- 
lution but noisy map and a featureless low resolution map. 
The weighting function is then related to the transform of 
this window function by 



w(9 g -e) = 



i 



(277)= 



T(k)J 2 (k.9)d 2 k 



(2.16) 



xi(0„-e) = 



(2.12) 



where J2 denotes the second order Bessel function. 

Finally, it is wor th re membering that the quantity a es- 
timated by equation (2.8) is the surface over-density in unit s 
of the critical surface density as defined by equation (2.6). 
Since Soxix depends on the geometry of the lensing configu- 
ration through the angular diameter-distance relationship, 
the mean redshift distribution of source galaxies is required 
before the absolute surface density can be obtained. In ad- 
dition, the KS technique is only sensitive to variations in 
surface density. This is because a uniform slab of material 
across the whole lens plane does not distort the images of 
galaxies lying behind. Thus, the mean surface density, S, is 
also unknown unless the region analysed is sufficiently large 
to encompass the whole of the lensing cluster so that the 
surface density near the edge of the region can be taken as 
the zero-point. 



and 



X2(0 g - 9) 



9x)(0 9 



(2.13) 



Note that both the ellipticities, e,, and the kernel, Xi, are P° _ 
lars, i. e. a rotation of the coordinate system through 180 de- 
grees leaves their components unchanged. This is perhaps 
most easily understood by visualizing each galaxy as an el- 
lipse which has 180 degree symmetry. 

The weighting function, W(9), is introduced to produce 
a smoothed estimate of the lens mass distribution. If no 



3 CREATING AND ANALYSING SIMULATED 
CCD IMAGES 

Our aim is to simulate B-band CCD images of lensed field 
galaxies. We assume that stars and the typically redder clus- 
ter galaxies have been identified and removed from the im- 
age. The CCD specification that we adopt is ~ 2000 by 
~ 2000 pixels each of size ~0.3 arcseconds on a side. These 
numbers were chosen to correspond to some observational 
data that we had obtained (see Wilson et al. 1995). The 
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magnitude limit we adopt is intended to correspond to what 
is currently achievable in one night on a 4-metre telescope. 

3.1 The Source Planes 

The distributions of galaxy size, ellipticity and redshift that 
we adopt are detailed below. These are intended to provide a 
realistic description of the mean distributions applicable to 
the majority of the galaxies which enter into the reconstruc- 
tion analysis. We simulate background galaxies with appar- 
ent magnitudes spanning the range 23.25 < me < 27.75, 
but the majority of the galaxies which enter our analysis 
have magnitudes close to the limit, m cu t = 25-26, which we 
impose in order to select a sample with well defined shapes 
and orientations (see Section 4.1). The inclusion of galaxies 
substantially fainter than m cut is required because the small 
proportion of these which happen to lie directly behind the 
cluster centre will undergo sufficient lensing amplification 
to fall subsequently within the detection limit. We adopt 
size, ellipticity and redshift distributions typical of those for 
galaxies with apparent magnitude ~ m cut and, for simplic- 
ity, we ignore variations in these distributions with apparent 
magnitude. We are less concerned here with modelling the 
genuine galaxy population which will, in all probability, be 
extremely complex, than with modelling a sensible but sim- 
ple population whose effects on the signal can more easily 
be followed. 

• Redshift and magnitude distributions 

As a reasonable redshift distribution for the simulated 
galaxies, we have adopted the me = 25 distribution pre- 
dicted by the analytic model of galaxy formation of Cole et 
al. (1994). As seen in Fig. 20 of that paper, the model has 
a median redshift of about z — 1, with a tail extending to 
z = 2.5. This model is broadly consistent with the observed 
redshift distributions of both bright and faint B- and K- se- 
lected samples. Since the critical density, 5 cr it, and therefore 
the bendi ng a ngle, j3, depend on the source redshift through 
equation (2JS), we discretely sample the redshift distribu- 
tion and produce a set of source planes spanning a range of 
redshifts. The net effect on the reconstruction is that 5 cr it 
is replaced by the mean value 



1 



■ p(z) dz 



(3.1) 



Scrit {z) 

where p(z) is the probability that a galaxy lies at redshift 
z. The distribution of apparent magnitudes we generate 
directly from the B-band source counts of Metcalfe et al. 
(1995). 

• Scale length distribution 

We assume that all the background galaxies are disks with 
exponential profiles, 



I(r) = 7 exp 



(3.2) 



This is a reasonable approximation since field galaxies are 
predominately disks. The scale length, A, we choose from a 
uniform distribution spanning the range 0.25 arcseconds to 
0.65 arcseconds, as suggested by observations (Tyson 1994). 
• Ellipticity distribution 

We take an empirical ellipticity distribution derived from 
a single 9.6 x 9.6 arcminute frame in 0.7-0.9 arcseconds see- 
ing on the 200 inch Hale telescope at Palomar (Brainerd, 



Blandford & Smail 1995). There are about 6000 galaxies 
catalogued in this frame and we sample their ellipticities 
at random. For a given magnitude and scale length, more 
elliptical galaxies have higher surface brightness. This is be- 
cause we assume conservation of intensity (i.e. no dimming 
by dust) and, since elliptical galaxies present a smaller cross- 
sectional area than face-on circular galaxies, they have a 
higher flux per unit area. Note that including a distribution 
of intrinsic ellipticity is equivalent to adding a noise term in 
the KS reconstruction procedure and this adds noise to the 
final surface over-density map. This is because intrinsic el- 
lipticity introduces scatter into the measured values of (ej). 
Any uncertainty in the (d) will translate into a correspond- 
ing uncer tain ty in the estimate of surface over-density via 
equation (2.8). 



We discuss the effects of varying these distributions in 
Section 5. 



3.2 The Lens 

In order to illustrate the effects of seeing and of non- 
linearities, we initially construct a simple spherically sym- 
metric lens with a Gaussian mass distribution. Later on, and 
for the main part of this investigation, we use a dark mat- 
ter cluster grown in a cosmological N-body simulation as a 
realistic complex lens. In each case we place the cluster at a 
redshift z = 0.18, which again corresponds to the observa- 
tions of Wilson et al. (1995). The N-body cluster, described 
in detail in Frenk et al. (1995), comes from a high resolution 
simulation of a cluster which was initially identified in a sim- 
ulation of a box 360 Mpc on a side of an Q = 1, Ho = 50 km 
s _1 Mpc _1 , cold dark matter universe (Davis et al. 1985). 
The initial conditions for the cluster were extracted from 
this large simulation and, after adding appropriate addi- 
tional high frequency noise, the cluster was simulated again 
with a P 3 M code using 262144 particles, this time in a box 
of size 45 Mpc. The spatial resolution in the simulation was 
35 kpc and the mass per particle 2.5 x 1O 1O M0. The cluster 
has a one- dimensional velocity dispersion of ~ 800 km s _1 . 

Our aim is now to calculate the bending angle on a 
grid of points corresponding to the centres of pixel s in the 
lens plane (Section 3.3). Using equations (2.4) and (2.5) we 
obtain the following relationship between the bending angle 
and the surface density 



V./3 



25 

5 C rit 



(3.3) 



It turns out that this equation is most easily solved in k- 
space, using Fast Fourier Transforms. In practice, we calcu- 
late the bending angle at a few points on a coarse grid and 
then use cubic splines to interpolate onto a finer grid. 



3.3 The Image Plane 

Next, we simulate the corresponding image plane. In the 
weak-lensing regime the bending angle a varies continuously 
and smoothly across the lens plane. Thus, to a good approx- 
imation, we can construct the image plane by mapping pixel 
by pixel the image plane from the source plane. The formula 
linking the positions of s our ce points and corresponding im- 
age points is equation (2^2). This formula allows 9 a to be 
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(1981) specifically for measuring properties of faint galax- 
ies. Galaxies have to satisfy certain criteria in order to be 
"detected". The user specifies an acceptable level of inten- 
sity and a minimum area. After detection FOCAS grows an 
isophote around each galaxy until it extends as far as the sky 
noise. The shape of the galaxy is then evaluated within this 
isophote. The values of intensity-weighted second moments 
from FOCAS are used to define the ellipticity components 
ei and ei that feed into equation (2.8) to yield the estimated 
surface over-density. 

The upper panel in Fig. 2 shows the surface over- 
density distribution of the N-body cluster evolved to red- 
shift z — 0.18, placed at the corresponding distance and 
smoothed with a Gaussian of am = 0.25 arcminutes. The 
lower panel shows, at the same resolution, the reconstructed 
surface over-density map obtained from our simulated CCD 
image in 1 arcsecond seeing. Although noise features are 
clearly visible in the reconstruction, the overall morphology 
is remarkably accurate, reproducing all the major features 
of the original cluster. We show these plots here as an illus- 
tration of the power of the method. In the next section we 
will investigate the accuracy of the reconstruction in more 
detail. 



4 THE RECONSTRUCTION METHOD IN 
PRACTICE 

In practice observations of gravitational lensing will not 
conform to the ideals assumed by the KS reconstruction 
method. The two most important limitations of real obser- 
vational data are seeing and noise. 



Figure 2. The upper panel shows the projected surface over- 
density of the N-body cluster which we use as the gravitational 
lens. The cluster is at a redshift of z = 0.18 and the surface over- 
density has been smoothed with a Gaussian of # sm = 0.25 arcmin- 
utes. At this redshift 1 arcminute corresponds to approximately 
0.24 Mpc. The lower panel shows the corresponding reconstruc- 
tion produced from a deep image in 1 arcsec seeing. In each panel 
the lowest contour corresponds to a = and the contour spacing 
is Act = 0.025. Note that a is the surface density in excess of the 
mean measur ed in units of the critical surface density for lensing 
(sec equation 2.6). For the redshift distribution assumed in this 



case S crit = 3.3 x 10 15 M Q /Mpc 2 

expressed uniquely in terms of Qi, but not Oi in terms of 8 S . 
For each image pixel we apply this formula and obtain the 
corresponding point in the source plane. We then assign an 
intensity to the image pixel by simple bilinear interpolation 
of the intensities in the nearest four source pixels. This pro- 
cedure results in a near perfect CCD image of the lensed 
galaxies. Finally, we add noise to the frame and then con- 
volve it with a Gaussian of width 8 see to simulate sky noise 
and atmospheric seeing. 



3.4 The Inversion 

We analyse the image built up in this way using FOCAS. 
FOCAS (Faint Object Classification and Analysis System) 
is a software reduction package developed by Jarvis et al. 



• Seeing 

Seeing is the distortion of images produced by scattering 
of light as it propagates through the Earth's atmosphere. 
Point sources become finite in extent and extended sources 
like galaxies undergo a corresponding blurring. The resulting 
effect is to make the galaxies appear more circular. This 
masks the true elongation of lensed galaxies, making the lens 
appear less strong and hence reducing the lensing surface 
over-density recovered by the KS method. 

• Noise 

Noise is any spurious signal introduced during the detec- 
tion process. There are various categories of noise e.g. pho- 
ton noise, background noise or detector noise. When observ- 
ing faint galaxies the most important source of noise is the 
sky background. The shapes of faint galaxies can be grossly 
distorted by noise. This can confuse the lensing analysis, so 
very faint galaxies need to be excluded (see Section 4.1). 

In addition the KS technique will also break down because 
of nonlinearity if the surface density of the lens is too high. 

• Nonlinearity 

The KS technique is applicable only to weak lensing sit- 
uations, when second order shear terms are negligibly small 
i.e. the bending angle varies only slowly, d/3/d6 <C 1. If the 
cluster surface over-density is large and varies rapidly this 
assumption is no longer valid and strong lensing techniques 
must be employed. 
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Figure 3. The ratio, /, of the true central surface over-density 
of the lens to the central surface over-density recovered by the 
KS technique. The upper panel shows the dependence of / on the 
central surface over-density of the lens for fixed seeing conditions 
of Osee — 1 arcsecond. The lower panel shows the variation of / 
with 8 BCC for a = 0.1. 

We illustrate the effect of seeing and non-linearity in 
Fig. 3. Here we use a simple spherically symmetric lens with 
a Gaussian mass profile. We vary both the seeing and the 
mass of the lens, but in each case we keep the noise added to 
the image frame at a very low level. The factor, /, plotted 
in both panels of Fig. 3 is the ratio of the true surface over- 
density, a, at the centre of the lens, to the corresponding 
value, a, recovered by the KS technique. The upper panel 
shows / as a function of the central surface over-density of 
the Gaussian lens, for a fixed seeing of 8 SCC — 1 arcsecond. 
We can see that / is constant up to about a = 0.1 and 
then increases with increasing surface over-density, implying 
that nonlinear effects are becoming important for surface 



over-densities greater than this value and the weak lensing 
approximation is beginning to fail. The lower panel shows / 
as a function of seeing. Here the central surface over-density 
of the Gaussian is kept fixed at a = 0.1 which, from the 
upper panel, is still in the regime where the weak lensing 
approximation appears valid. We can see that / increases 
rapidly as the seeing worsens. 

It is not easy to correct for the systematic error caused 
by non-linearity. Thus, if the surface over-density at the cen- 
tre of massive clusters is to be accurately estimated, alterna- 
tive methods must be employed which take account of non- 
linearity (C. Seitz & Schneider 1995a; Kaiser 1995). Note, 
however, that the systematic error is only of order 25% at 
the centre of a cluster of central surface over-density 0.3. 
(Note also that the reason why / does not tend to unity at 
low values of a, in the upper panel of Fig. 3, is the 1 arc- 
second seeing, not residual non-linearity.) Elsewhere in the 
cluster the systematic error will be smaller. 

The systematic error due to the blurring effect of seeing 
is potentially much larger. If we were able to tabulate the 
ratio / for all observing conditions, then this table could 
be used to find a compensation factor to correct the surface 
over-density estimates returned by the KS technique. How- 
ever, this is not practical since the effect of seeing depends 
not only on the value of # see , but also on intrinsic properties 
of the galaxy images used in the reconstruction. For exam- 
ple, the degradation due to seeing increases as the angular 
size of the galaxy images used decreases. To circumvent this 
problem we outline a calibration procedure in Section 4.2 
which can be used to estimate the required compensation 
factor, /, for any given observational dataset. 

4.1 Defining the Galaxy Sample 

In our simulations we have assumed that stars and clus- 
ter galaxies have been removed from the CCD image. In 
practice, since cluster galaxies are mostly E/SO's and are all 
at approximately the same redshift, they have very similar 
colours and, provided two colour information is available, 
they are relatively easy to identify. On a colour-magnitude 
diagram of all objects within the frame, the cluster galaxies 
will fall on a (nearly) horizontal line (see e.g. Smail 1993) 
and can be excluded from any subsequent lensing analysis. 

The value of the signal-to-noise ratio in our simulations 
has been chosen to mimic detections of galaxies down to 
m,B — 26.5. Although all galaxies down to this magnitude 
limit are detected in the simulated CCD frame, the faint 
galaxy shapes are badly contaminated by noise. Thus, it is 
necessary to make a cut at a brighter magnitude in order to 
exclude these faint galaxies from the analysis. We find that a 
useful guide to selecting this magnitude cut comes from look- 
ing at the ellipticity distribution of the galaxy images as a 
function of apparent magnitud e. It is t o be expected that the 
intrinsic distribution of e = yef + ef will be a slowly vary- 
ing function of apparent magnitude. It is, in fact, assumed 
to be constant in our simulations. Thus, a sudden change in 
the shape of this distribution at faint magnitudes can most 
likely be attributed to the onset of noise in the image cor- 
rupting the shapes of the faint galaxies. As a quantitative 
comparison, we divide the data into half magnitude bins 
and compare the ellipticity distribution in each bin in turn 
with a representative bright sample using the Kolmogorov- 
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Figure 4. The probability that galaxies in any given half- 
magnitude subsample are drawn from the same distribution of 
cllipticities as that of the 23.75 to 24.25 subsample. The proba- 
bilities are calculated by means of a Kolmogorov-Smirnoff com- 
parison test. Populations fainter than m=25.5 are inconsistent 
with the brighter populations. 



Smirnoff comparison test. As can be seen from Fig. 4, the 
probability that the two distributions are drawn from the 
same population plummets to virtually zero at magnitudes 
fainter than m=25.5. We have found this transition to be a 
good indication of where to place the magnitude cut used to 
define the sample of galaxy images to be fed through the KS 
reconstruction technique. The appropriate value of the cut 
will depend, of course, on the specific observational set-up. 

4.2 Calibration of mass estimator 

In this section we describe how to calibrate a CCD frame 
for use in the KS reconstruction method. Specifically, we 
show how to compute a compensation factor, /, that corrects 
for the bias in the surface over-density estimates returned 
by the KS method. Briefly, the procedure involves shearing 
the galaxy images by a known amount, adding seeing and 
measuring the resultant shear. The compensation factor, /, 
is then the ratio of the input shear to the measured shear. 

If one takes the image frame, multiplies the x-coordinate 
of each pixel by a factor 1 + and rebins, then all the galaxy 
images will be sheared in the x-direction. If e is small and 
the initial distribution of ellipticities is not too broa d then it 
is easy to show from the definitions ( ) and ( ^.10|) that the 
ellipticity component e% of each galaxy is unchanged while 
the ei component is on average increased by e. If the galaxy 
images are then blurred by seeing one will find that the 



Figure 5. The compensation factor, /, versus seeing in realistic 
signal-to-noise conditions. The line shows the values obtained for 
e cut = 0.5 



measured change in the shear will be somewhat less than e. 
Since, according to equation (2.8), the surface over-density 
at any given point is proportional to the measured elliptici- 
ties, the ratio of e to the mean change in ei, (Aei), is in fact 
the factor / required to correct the surface over-densities, 



/-v = 77 f T (4.1) 
a (Aei) 

This procedure is complicated by two factors. First, the 
value of / depends on the sizes of the galaxy images (A given 
value of seeing will produce a much greater circularising ef- 
fect on small galaxies than on large galaxies). Hence one 
would underestimate f if the shearing process were applied 
directly to the enlarged blurred observed images. Second, 
the initial distribution of ei can be quite broad with some 
galaxies having values of ei approaching unity prior to ad- 
dition of any further shear. Since e\ is constrained to be 
less than unity, these high values of ei cannot be increased 
further. It is therefore necessary both to deconvolve the im- 
age prior to applying the shear and to limit the analysis to 
galaxies whose original ellipticity is less than some value, 
e cu t • 

In summary, our calibration procedure consists of the 
following steps: 

(i) Deconvolve the CCD image using the Point Spread 
Function measured from one or more stars on the frame. 
Note that since no analysis is to be made using the decon- 
volved images, it is not necessary to use sophisticated noise 
suppressing deconvolution algorithms. Neither is it required 
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to model the PSF in great detail. A Gaussian fitted to the 
measured PSF is probably adequate. 

(ii) Stretch the galaxy images along the x-axis by a known 
factor, 1 + e, and rebin. A value of e « 0.1 is appropriate as 
this is typical of the values produced by weak lensing. 

(iii) Reconvolve the stretched image with the same PSF. 

(iv) Run reduction software on both the original image 
and this new stretched image and compute the ellipticity 
components ei and ei for each galaxy. 

(v) Select the galaxies with measured values of e < e cut in 
the original frame and, for these, compute the mean change 
in ei, (Aei), between the original and stretched frames. De- 
fine the compensation factor / = e/(Aei). 

(vi) Estimate th e le ns surface over-density using the KS 
method, equation (2.8), with the galaxies selected using the 
same cut in e as above. Finally, multiply the resulting surface 
over-densities, a, by the factor / to yield corrected estimates. 

Fig. 5 shows estimates of / obtained by this proce- 
dure from simulated CCD frames constructed with realistic 
signal-to-noise ratios and galaxy populations. Each point 
is the average obtained from 10 simulated CCD images. 
The curve shows the compensation factor / for the choice 
e C ut = 0.5. As expected, we see that / increases rapidly as 
seeing worsens and that even for good seeing conditions it is 
significantly different from unity. We use values of / calcu- 
lated in this way in the next section and show that they give 
remarkably good estimates of the true surface over-density. 



5 EXAMPLES AND DISCUSSION OF 
RESULTS 

We now examine a series of test cases where we explore 
the success and reliability of the calibrated reconstruction 
technique for a range of observational conditions and sample 
selections. 
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Figure 6. For simulated observations taken in conditions of 
Sscc = 0.6 arcseconds seeing the panels show the reconstructed 
map of the lens surface over-density, a scatter plot of the com- 
pensated estimated versus true surface over-density and the esti- 
mated shear pattern across the face of the lens, respectively. The 
value of the compensation factor, /, estimated from this data is 
also shown on the middle panel. The details of the assumed back- 
ground galaxy properties are detailed in Section 3.1. A magnitude 
cut, r?icut = 25.5, suggested by the Kolmogorov-Smirnoff test of 
Section 4.1, and an ellipticity cut of e cu t = 0.5 were used to define 
the galaxy sample that was analysed. 
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Figure 7. As Fig. 6 but for seeing of I 



1.0 arcseconds. Figure 8. As Figs. 6 and 7, but for seeing of 6 S , 

onds. 



1.6 arcsec- 



5.1 The Effect of Seeing 

Figures 6,7 and 8 show the results of the calibrated recon- 
struction technique for three different values of the seeing, 
# S ec = 0.6, 1.0 and 1.6 arcseconds. The upper panel in each 
figure shows the reconstruction of the cluster surface over- 
density map. This should be compared to the true lens sur- 



face over-density displayed in Fig. 2. The central panel is 
a scatter plot of the compensated estimated surface over- 
density versus the true surface over-density measured on a 
62 x 62 grid covering the region shown in the upper panel. 
The compensation factor / estimated as described in Sec- 
tion 4.2 is shown in the upper left of the panel. The lower 
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Figure 9. The reconstructed surface over-density map and scat- 
ter plot for 9 Bee = 1.0 arcseconds (as in Fig. 7), but for galaxies 
with scale lengths 20% larger than assumed in Fig. 7. 



Figure 10. As Fig. 7, but for a narrower distribution of intrinsic 
cllipticitics for the background galaxies. The mean of b/a is now 
0.79 compared with 0.71 in Fig 7. 



panel depicts the mean shear of the galaxies in 10 x 10 bins, 
again covering the same area as the map in the upper panel. 
The length of each line is proportional to the mean elliptic- 
ity, e = y (ei) 2 + (e2} 2 , in each cell and the orientations of 
the lines indicate the direction of the shear. 

The first point of note is that the complex morphol- 
ogy of the cluster mass distribution is recovered quite well 
in all three cases, with only a gradual degradation of the 
reconstruction as the seeing becomes progressively worse. 
Second, the value of the compensation factor, /, is a strong 
function of the seeing, varying from / = 1.45 for 8 SCC — 0.6 
to / = 2.21 for f) scc = 1.6. In spite of this, the compensated 
surface over-density estimates, fa, are in good agreement 
with the true values, a, i.e. the points in the central pan- 
els all scatter around the line fa = a , wit h no significant 
bias. Finally, we note that equation ( 2.15 ) appears to be 
a good approximation to the variance in the surface over- 
density estimator; from 9 simulations in 1 arcsec seeing we 
find a scatter in a of 0.052 which compares well with the 
theoretical value of 0.045. 



5.2 Variations in the Properties of the 
Background Galaxies 

In our simulations so far we have assumed specific distribu- 
tions of the sizes, shapes and redshifts of the background 
galaxies. These are realistic examples but it is nevertheless 
important to investigate how our results change when they 
are varied. In Figs. 9 to 12 we explore the effect on the re- 
constructed surface over-density maps and scatter plots of 
varying each of these distributions in turn. In all cases we 
employ 1 arcsecond seeing and therefore compare our results 
to Fig. 7. The factors / are now calculated from a single 
frame rather than from the mean of ten frames as before, so 
they will be somewhat less reliable. 

In Fig. 9 we experiment with the sizes of the galaxies 
by assuming they are 20 percent larger than before. In this 
case we might expect the scatter in fa to remain unchanged, 
but / to be reduced since the ratio of galaxy size to seeing 
is increased. Indeed, the scatter in Fig. 9 is similar to that 
in Fig. 7 and the value of f is reduced. The scatter does not 
change because the ellipticity distribution and the noise, the 
major contributors to the uncertainty, are unchanged. 

In Fig. 10 we explore the effect of varying the distribu- 
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Figure 11. As Fig. 7, but with the background galaxies all at 
redshift z = 2. 



Figure 12. As Fig. 
galaxies down to m cu 



but for a deeper CCD frame in which 
= 26.5 are retained. 



tion of intrinsic galaxy ellipticities. We bias the distribution 
slightly towards less elliptical galaxies, so that the mean ax- 
ial ratio is b/a = 0.79 rather than b/a = 0.71 as before. 
We see from the figures that the resultin g eff ect is to reduce 
the scatter in accordance with equation (2.15). Although the 
value of f shown here is slightly diminished, the mean value 
of f from 5 frames is little changed. 

In Fig. 11 we use a deeper distribution of galaxy red- 
shifts, namely placing all the galaxies at z = 2. This change 
reduces the value of 5 cr it (see equation 2.6) and increases the 



values of a and a proportionally at all grid points. This is 
the reason for the change of scale on the axes compared with 
the previous figures. Since the values of a and a increase in 
the same ratio, the mean value of / is unchanged. 

In our final plot, Fig. 12, we increase the signal-to- 
noise ratio so that galaxies can now be detected one magni- 
tude fainter than before. The higher galaxy number density 
great ly re duces the scatter in fa, as expected from equa- 
tion ( j2.15| ). The compensation factor is also reduced some- 
what because the galaxy shapes are now less distorted by 
noise than before. 



6 CONCLUSIONS 

We have performed a series of controlled experiments to as- 
sess the reliability of the technique proposed by Kaiser & 
Squires (1993) to reconstruct the surface mass over-density 
of galaxy clusters from observations of weak gravitational 
lensing. In particular, we have tested the KS method on a 
realistic cluster mass distribution, typical of those expected 
in an Q. — 1 universe. By simulating data from standard ob- 
serving conditions, we have investigated the effects of seeing 
and signal-to-noise on the reconstructed dark matter maps 
and we have explored how the results vary with different 
assumptions for the distributions of intrinsic galaxy shapes 
and redshifts. Our main conclusions are as follows: 

1. With a careful analysis of data obtained in standard ob- 
serving conditions, the KS method provides a remarkably 
faithful reconstruction of the morphology of a complex clus- 
ter, reproducing the richness of structure expected in an 
fl = 1 universe. 

2. Our simulations show that the weak lensing assumption 
on which the KS technique is based begins to break down 
when the mass surface over-density exceeds 10% of the criti- 
cal surface density. However, even when it equals 30% of the 
critical value, the KS method underestimates the surface 
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over-density by only ~ 25%. The noise in the recon struc ted 



maps agrees well with a simple estimate (equation 2.15) of 



the uncertainties due to Poisson noise and to the scatter in 
the intrinsic ellipticities of the lensed background galaxies. 

3. The simple calibration procedure for CCD images which 
we have designed and tested, efficiently corrects for the ef- 
fects of atmospheric seeing on the reconstructed mass sur- 
face over-density maps. This procedure is straightforward to 
apply to CCD data and yields a multiplicative "compensa- 
tion factor", /, which allows the true surface over-density in 
each pixel (in units of the critical density) to be recovered 
from the reconstructed map. The absolute value of the lens 
mass cannot be derived by this method unless the critical 
density, which depends on the redshift distribution of the 
lensed galaxies, is known. However, the dependence on red- 
shift is fairly weak (see e.g. Fig 5 of Blandford & Kochanek, 
1987). 

4- Our method for calculating the compensation factor is 
quite robust. The value of / is primarily determined by the 
seeing, but the depth of the CCD image and the intrinsic 
properties of the lensed galaxies also affect it. Useful results 
can be obtained even with data acquired in seeing as large 
as 1.6 arcseconds, although the technique clearly works best 
with sub-arcsecond seeing. Our simulations indicate that it 
should work well when applied to data obtained in a wide 
range of observing conditions, allowing useful mass recon- 
structions to be made even when the correction factor is as 
large as / = 2. 
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